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ABSTRACT 

Aims. The knowledge of accurate stellar parameters is a key stone in several fields of stellar astrophysics, such as asteroseismology 
and stellar evolution. Although the parameters can be derived both via spectroscopy and multicolor photometry, the obtained results 
are sometimes affected by systematic uncertainties. In this paper we present a self-consistent spectral analysis of the solar- type star 
HD 49933, which is a primary target for the CoRoT satellite. 

Methods. We used high-resolution and high signal-to-noise ratio spectra to carry out a consistent parameter estimation and abundance 
analysis of HD 49933. The LLmodels code was employed for model atmosphere calculations, while SYNTH3 and WIDTH9 codes 
were used for line profile calculation and LTE abundance analysis. 

Results. In this paper we provide a detailed description of the methodology adopted to derive the fundamental parameters and the 
abundances. Although the obtained parameters differ from the ones previously derived by other authors, we show that only the set 
obtained in this work is able to fit the observed spectrum accurately. In particular, the surface gravity was adjusted to fit pressure- 
sensitive spectral features. 

Conclusions. We confirm the importance of a consistent analysis of relevant spectroscopic features, application of advanced model 
atmospheres, and the use of up-to-date atomic line data for the determination of stellar parameters. These results are crucial for further 
studies, e.g. detailed theoretical modelling of the observed pulsation frequencies. 
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1. Introduction 

In the last ten years several space missions were launched to ob- 
tain high quality photometric data (e.g. WIRE, MOST, CoRoT) 
and several others are going to be launched in the near future 
(e.g. KEPLER, BRITE). One of the main goals of these mis- 
sions is to provide astronomers with high precision photometric 
data that will allow a better understanding of stellar pulsation - 
the only approach that can improve our knowledge about stellar 
interiors. 

The modelling of pulsational signals requires the knowl- 
edge of stellar parameters and primarily accurate values of the 
effective temperature (T e g) and metallicity (Z). The determi- 
nation of fundamental parameters can be performed by differ- 
ent methods (some examples for stars from B- to G-type are : 
iFuhrmann etal]|1997t iPrzvbilla et afll2006t iFossati et alj|2009l) 
that not always lead to consistent results. Thus it is important to 
choose a methodology that allows to constrain the parameters of 
the star from the available observables (usually photometry and 
spectroscopy) in the most robust and reliable way. 

After Moss eTet al.l (|2005) discovered the presence of solar- 
like oscillations in HD 49933 using RV measurements from high 
resolution time-series spectra, this star was included as a pri- 
mary CoRoT target for the photometric observation of solar- 
like oscillations. It was one of the first main sequence solar-type 
star in which solar-like oscillations were clearly detected from 



space photometry, t hanks to the high qualit y data provided by 



l qu; 

the CoRoT satellite dAppourchaux et al.|[2008b . 

Fundamental parameters and abundances of HD 49933 were 
obtained previously by different authors with different meth- 
ods. The published fundamental parameters of HD 49933 dis- 
play some scatter both in T e s a nd log g (e.g. lEdvardsson et al.l 
1 19931: iBlackwell &Lynas-Gravl 11998 billon & Magainl 12006: 
ICenarro et al.l 120071 : iBruntt etaU 120081) . As a consequence, 
the published abundances show some scatter. At the same 
time a need for a higher precision in stellar parameters for 
HD 49933 was e mphasised by recent asteroseismic mod- 
elling of this object ( Appourcha ux et ail 120081 : Kallinger et al.l 
2009 ), which employed new pulsation analysis techniques (e.g. 
bruberbauer et alj|2008t) . 

The main goal of the present work is to perform a consis- 
tent atmospheric and abundance analysis of HD 49933 which 
fits all photometric and spectroscopic data and to investigate to 
what degree the derived fundamental parameters depend on the 
applied methods. 



* The CoRoT space mission was developed and is operated by the 
French space agency CNES, with participation of ESA's RSSD and 
Science Programmes, Austria, Belgium, Brazil, Germany, and Spain 



2. Observations and spectral reduction 

Ten spectra of HD 49933 were obtained between February 1 1th 
and February 13th 2006 with the cross-dispersed echelle spectro- 
graph HARPS (spectral resolution R ~ 115 000) attached to the 
3.6-m ESO La Silla telescope. The spectra were reduced with 
the online pipelin^B 



1 http : //www . eso . org/sci/f acilities/lasilla/instruments/harps/t( 
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Each spectrum was obtained with an exposure time between 
243 and 300 seconds. We retrieved the spectra from the ESO 
archive and co-added them to obtain a single spectrum with a 
signal-to-noise ratio (SNR) per pixel of about 500, calculated 
at ~5000 A. The final spectrum, normalised by fitting a low order 
spline to carefully selected continuum points, covers the wave- 
length range 3780-6910A. The spectrum has a gap between 
5300 A and 5330 A, since one echelle order is lost in the gap 
between the two chips of the CCD mosaic detector. 

Figure|4](Online material) shows (from top to bottom) a sam- 
ple of the first spectrum obtained on February 1 1th, the last ob- 
tained on February 13th and the final co-added spectrum, around 
the strong Fe n line at 5018 A. 

3. The model atmosphere 

To compute model atmospheres of HD 499 33 we employed th e 
LLmodels stellar model atmosphere code (Shul yak et al.l l2004). 
For all the calculations Local Thermodynamical Equilibrium 
(LTE) and plane-par allel geometry were assumed. We used 
the VALD database dpTskunov et all fl995l iKupka et all 1 19991: 
Ryabchik ova et al. l l 19991) as a source of atomic line parameters 
for opacity calculations. For a given model atmosphere we per- 
formed a line selection procedure which allowed to chose lines 
that significantly contribute to the opacity for a given set of 
model parameters, adopting the selection threshold £ v /a v ^ 1%, 
where a v and l v are the continuum and line absorption coeffi- 
cients at a giv en frequency v. Convection wa s implemented ac- 
cording to the ICanuto & Mazzitelhl (Il991albl) model of convec- 
tion (see iHeiter et aLll2002i for more details). 



4. Fundamental parameters and abundance 
analysis 

iGillon & Magarnl ([2006) 

used Stromgren indices to deter- 
mine the atmospheric para meter s of H D 49933. By using 
the TEMPLOGG code (iRogersI 1 19951) they found r eff = 
6543+200 K, logg = 4.24±0.20and [Fe/H] = -0.38 + 0.20 dex. 
Although the uncertainties are quite large, we used this as our 
starting point and in an iterative process to gradually improve 
the parameters using different spectroscopic indicators as we 
will describe in detail below. In our analysis, every time any 
of r e ff, logg, i/ m ic or abundances changed during the iteration 
process, we recalculated a new model with the implementation 
of the last measured quantities. The same concerns newly de- 
rived abundances: while the results of the abundance analysis 
depend upon the assumed model atmosphere, the atmospheric 
temperature-pressure structure itself depends upon the adopted 
abundances, so that we recalculated the model atmosphere every 
time abundances were changed, even if the other model parame- 
ters were fixed. Such a procedure ensures the model structure to 
be consistent with the assumed abundances. 

We performed the T e s determi nation by fitting sy nthetic line 
profiles, calculated with SYNTH3 (lKochukhovl2007l) . to the ob- 
served profiles of three hydrogen lines: Ho;, H/3 and Hy. In the 
temperature range expected for HD 49933 hydrogen lines are ex- 
tremely sensitive to temperature variations and very little to log g 
variations, therefore being good temperature indicators. The 
obtained with this procedure is T e ff = 6500 + 50 K. 

Figure[T]shows the comparison between the observed H/3 line 
profile and the synthetic profiles calculated with the adopted stel- 
lar parameters, as well as the synthetic p rofile obtained wit h the 
higher r eff = 6780 + 130K published by|B^untLeLaL|(l2008j). 



observation 

--■ Teff=6500K 
Teff=6780K 




Wavelength (A) 

Fig. 1. Comparison between the observed H/3 line profile (solid 
line) and synthetic profiles calculated with the final adopted 
r eff = 6500 K (da shed line) and with the r eff =6780K from 
iBruntt et~ai] (Hp08) (solid thin line). The dashed line agrees al- 
most perfectly with the observed spectrum. 



Another spectroscopic indicator for T e g is given by the anal- 
ysis of metallic lines. In particular, r e ff is determined eliminat- 
ing the correlation between line abundance and the line excita- 
tion potential Ofexdt) for a given ion/element. This procedure can 
lead to erroneous parameters, so we decided not to take into ac- 
count this indicator in our analysis. This important point will be 
extensively discussed in Sect. [5] 

The surface gravity was derived from two independent meth- 
ods based on line profile fitting of Mgi lines with developed 
wings and ionisation b alance for several ele ments. The first 
method is described in Fuhr mann et a i] (ll997h and is based on 
the fact that the wings of the Mgi lines at AA 5167, 5172 and 
5183 A are very sensitive to logg variations. In practice we de- 
rived first the Mg abundance from other Mg i lines without de- 
veloped wings, such as AA 4571, 4730, 5528 and 57 11 A, and 
then we fit the wings of the other three lines by tuning the 
log g value. To apply this method, it is necessary to have very 
accurate log gf values and Van der Waals (log ywaais) damp- 
ing constants for all the lines. The log gf values of the Mg 
lines given in the VALD database are of rather high quality 
and came from the laboratory measurements of the lifetimes 
dAnderson et alj[l967l) . However, recently new laboratory mea- 
surements of the oscillator strengths for Mgi lines at AA 5167, 
5172 and 518 3 A with an ac c uracy crlog gf- +0.04dex were 
published by lAldenius et alj d2007l). Van der Waals damping 
constants in VALD are calculated bv lBarklem & O'Maral (2000) 
but they appear to be slightly higher tha n needed to fit the 
solar li nes. While logywaais = -7.27 in iBarklem & O' Mara 
d2000h , lFuhrmann et alj d 19971) derived l ogywaais - -7.42 from 
the fitting of the solar lines using the lAnderson et al. I d!967h 
oscillator strengths. For our analysis we employed different 
combinations of oscillator strengt hs and damping cons tants. 
Adopting oscillator stre ngths from | Aldenius et alj d2007l) and 
damping constants from Barklem & O'Mara (2000) or oscilla- 
tor st rengths fro m Anderson et al.l d 1967b and damping constants 
from Fuhrman n et al.l d 1997b we derived a logg value of 3.85 , 
while adopting oscillator strengths from Aldenius et"aT] d2007l) 
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Fig. 2. Comparison between the observed Mg i line profiles used 
to determine log g (thick line) and synthetic profiles for the Sun, 
Procyon and HD 49933, from top to bottom. The synthetic spec- 
trum shown for HD 49933 was calculated with the final adopted 
logg= 4.00 (thin line). For all three stars we adopted the same 
combination of log gf and Van der Waals dam ping constants: os- 
cillator stre ngths from | Aldenius et all (l2007h and damping con- 
stants from iFuhrmann etall (1 19971) . See text for details on the 
observed spectra and adopted model atmospheres for the Sun 
and Procyon. The spectra of the Sun and Procyon were shifted 
upwards by 1 and 0.5 respectively. 



and damping constants from Fuhr mann et al.l d 1997b we obtained 
a logg value of 4.00. We obtained a formal logg value of 3.93 
+ 0.07 where the given error bar on logg mainly depends on 
the uncertainty of the damping constant (ywaais)- Comparison of 
the synthetic colors with the observed ones (see Sect. |5.21 > favors 
a logg value of 4.00 + 0.15. Figure [2] shows a comparison be- 
tween the synthetic and observed Mgi line profile for the 5172 
and 5183 A M gi lines. The same procedure applied to the solar 
flux spectrum (iKurucz et alj [l984) resulted in log g-4.44 (oscil- 
lator strengths from Aldenius et al. 2007 and damping constants 
from Fuhrmann et al. 1997) an d log #=4.30 for the c ombination 
of the oscillator strengths from Aldenius et al. (2007) and damp- 
ing constants from lBarklem & (TM ara (2000). The comparison 
of the observed and the synthetic Mgi line profile for the 5172 
and 5183 A Mgi lines for the Sun is shown in Fig. [2] 

For surface gravity determination one often uses ionisation 
equilibrium, but this method is extremely sensitive to the non- 
LTE effects present for each ion/element, while Mg lines with 
developed wings (less sen sitive to non-LTE effe cts) are more 
suitable as logg indicators dZhao & Gehren 2000). For this rea- 
son we decided to keep the Mg lines as our primary log g in- 
dicator checking afterwards the obtained result with the ionisa- 
tion equilibrium. The log g value we obtained (log g=4.00+0. 15) 
is l ower than that given by photometry and by o ther authors, 
_ Gillon & Magainl 12006): log #=4.26+0.08 and lBruntt etaT] 
1) : logff=4.24±0.13 . The difference be t ween o ur result and, 
for ex ample, the one of [Gillon & Magain (2006); B runtt et ail 
(2008) is clearly connected with the difference in the derived ef- 
fective temperature, while there is a reasonable agreement for 
logg, but the error bars are unfortunately quite large. This im- 
portant point will be discussed in Sect. [5] 

Since HD 49933 does not have an effective temperature high 
enough to show He lines, we are not able to measure the at- 



mospheric He abundance. Diffusion calculations for low-mass 
metal poor stars s how that helium should b e depleted in the at- 
mospheric layers dProffitt & Michaudl 19911) . therefore we tested 
the effect of a helium underabundance on the obtained parame- 
ters. We rederived T e s and log g respectively with hydrogen lines 
and Mg i line profiles assuming a helium abundance of -2.0 dex 
and of -4.0 dex. In both cases we obtained no visible change 
of the hydrogen line profiles indicating no effect on the T e g de- 
termination, while logg should be increased of about 0.15 with 
both helium underabundances, to be able to fit again the Mg i line 
profiles adopted to derive log g. This value is within the adopted 
uncertainty for log g. 

Our main source for the atomic parameters of spectral lines 
is the VALD database. LTE abundance analysis was based on 
equiv alent widths, analysed with a modifie d version (Tsvmbal 
1996) of the WIDTH9 code dKuruczlll993ah . Although 301 Fei 
lines were measured, only 158 lines with accurate experimental 
oscillator strengths were chosen for abundance determination to 
achieve the highest accuracy. First, we reject the lines with the- 
oretical oscillator strength s. The rest of the fi nes were checked 
in the solar flux spectrum dKurucz et alJI 19841) . observed with a 
resolving power R 340000 at wavelengths between 4000 A 
and 4700 A and with R = 520000 at longer wavelengths. The 
final choice was made for the lines that did not require an os- 
cillator strength correction greater than 0. 1 dex to fit the line 
cores of the solar spectrum. We also tried to have a set of Fe i 
lines uniformly distributed over the range of equivalent widths 
and excitation potentials. As for the lines of other elements/ions 
we used nearly all unblended spectral lines with accurate atomic 
parameters available in the wavelength range 3850-6880 A, ex- 
cept lines in spectral regions where the continuum normalisa- 
tion was too uncertain. In case of blended lines, for lines sub- 
jected to hyperfine splitting (hfs) or for lines situated in the wings 
of the hydrogen lines we derived the line abundance perform- 
ing synthetic spectrum calculations with the SYNTH3 code. 
The hfs constants for abun d ance c alculations were taken from 
iBlackwell- Whitehead et al.l d2005l) for Mn i lines, from |Biehll 
d 19761) for Ban A 4554 A fine and from lLawler et alj d2001l) for 
Eun fines, hfs effects are negligible for the two Cui lines we 
used in the abundance analysis. 

The projected rotational velocity and macroturbulence 
(i/ macro ) were determined fitting synthetic spectra of several care- 
fully selected lines in the observed spectrum. We obtained v sin i 
= 10 ± 0.5kms~' and i/ macro = 5.2 + 0.5kms~'. The value de- 
rived for i/ maC ro is in agreement with th at expected according to 
the relation T e ff - i» macro published by IValenti & Fished {2005) 
and obtained on a sample of more than a thousand stars. 

Moss er et al. I d2005l) detected a line profile distortion variable 
with time with an amplitude of less than 500 m s _1 , probably due 
to granulation. This distortion is too small to be able to affect 
the parameter and abundance determination since the parameters 
were mainly derived with hydrogen and Mg i line profiles, while 
abundances used equivalent widths. The distortion of the line 
profiles could have, at most, increased the microturbulence ve- 
locity. It is known that Am stars show also distorted line profiles 
that lead to an increase in v m i c of 1.5-2. 0km s -1 , but for Am s tars 
the distortion is of the order of 3-4 km s" 1 dLandstreetlll998l) . If 
granulation played a role in increasing the derived value of v^c, 
it is probably less than 0.1 kms~', which is below our detection 
limit. 

The microturbulence was determined by minimising the cor- 
relation between equivalent width and abundance for several 
ions. We mainly used Fei lines since this is the ion that pro- 
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Table 1. LTE atmospheric abundances HD 49933 with error estimates based on the interna l scatter from t he n umber of analysed 
lines, n. Fourth and fifth columns give HD 49933 abundances relative to the solar values from lAsplund et ail (120051) (AGS) and from 
iGrevesse & Sauva"fl(ll998l) (GS), respectively. The sixth column gives abundances der ived by Bruntt et al. (2008) relative to GS. The 
last column gives the abundances of the solar atmosphere from lAsplund et ail d2005l) . 
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vides the largest number of measured lines within a wide range 
in equivalent width, but also the correlations obtained with Ti i, 
Ti n, Cri, Fe n, and Ni i lines were taken into account. The range 
of the microturbulent velocities goes from 1.4kms _I (Tii) to 
1.9kms _I (Tin) with an average of i; m i c =1.60 + 0.18kms _1 . 
The microturbulent velocity derived from Fe i lines is 1 .5 km s . 
This value, with an uncertainty of 0.2kms _I , was adopted for 
the final analysis as the best representation for all lines of neu- 
tral elements. Figure[5](Online material) displays the correlation 
between the line abundance of Fe i and the measured equivalent 
widths. 

The final abundances are given in TableQ] We also computed 
the abundance differ ence between HD 499 33 and the solar atmo- 
sphere as derived bylAsplund et alj d2005) (4th column) and by 
IGrevesse & Sauval d 19981) (5th column). The results of the abun- 
dance analysis of HD 49933 taken from lBruntt et alj d2008l) are 
given in the 6th column of Table[T]for comparison. 



The stellar metallicity (Z) is defined as follows: 



Zctoi- 



•2 IQLogiNJNu,) 



Z^mJOW^)' 



(1) 



where a is the atomic number of an element with atomic mass 
m a . Making use of the abundances obtained from the performed 
analysis we derived a metallicity of Z - 0.008 + 0.0 02 dex, 
adopting the solar abundances by lAsplund et al.1 d2005l) for all 
the elements that were not analysed. However, if we assume 
an underabundance of -0.5 dex for all not analysed elements 
the resulting Z value remains practically unchanged. Scaling the 
solar abundances of all elements by -0.5 dex gives the metal- 
licity Z = 0.006 dex. This substantial difference in Z illustrates 
that it is important to have an accurate determination of C and 
O, which make a large contribution to the determination of Z. 
Additionally, a He underabundance of - 1 and -3 dex relative to 
Sun does not change the Z value. 
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4. 1 . Abundance uncertainties 

The abundance uncertainties shown in Table Q] are the standard 
deviation of the mean abundance from the individual line abun- 
dances and do not represent the real error bars associated with 
each derived element abundance. 

Based on realistic errors on the equivalent width measure- 
ments we find that the uncertai nty on the abundan ces of a given 
line is only about 0.01 dex (see Fossati et al. 2009, for more de- 
tails). In the case of ions with a sufficiently high number of lines 
we assume that the internal scatter for each ion includes the un- 
certainties due to equivalent width measurement and continuum 
normalisation. 

Plotting the abundance scatter as a function of the number 
of lines it is possible to conclude about the mean scatter that 
can be applied in cases when only one line of a certain ion is 
measured. In these cases it is reasonable to assume an internal 
error of 0.10 dex. 

The internal scatter is just a part of the total abundance error 
bar since the uncertainties on the fundamental parameters are 
also playing an important role. Table [3] (Online material) shows 
the variation in abundance resulting from the changing of each 
parameter (T e g, \ogg and v m \ c ) by +lcr and the final adopted 
uncertainty for each ion following the standard error analysis. 

The main source of the uncertainty is the effective tempera- 
ture, while the uncertainties due to log g and in particular to v^c 
are almost negligible. 

Assuming that the contributions to the uncertainty on the 
abundances are independent, we derive the final uncertainty us- 
ing error propagation. The result is given in column seven of 
Table [3] (available as Online material). 



5. Discussion 

5.1. Comparison with previous determinations 

Our abundance calculation with the use of different atmospheric 
models for HD 49933 proposed in literature shows that the accu- 
racy (rms) of the abundance determination is the same for prac- 
tically all of them and any difference in absolute value is caused 
by the differences in the model parameters. Therefore we shall 
compare and discuss here the determination of the atmospheric 
parameters. 

The publi shed effective temperatu re values range 
from 6300 K dHartmann & Gehrenl Il988l) through 6484 K 
Blac kwell & Lvnas-Gravl I1998I) (IR flux method), 6595 K 
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(Edvardsson et al. 1993) (Stromgren photometry ) to 6 735- 
6780 K (iGillon & Magainll2006l IBruntt et~ai]l2004 [2008b . The 
last three cited works employed spectrum synthesis iterative 
methods: APASS, VWA. For more details o n previously derive d 
atmospheric parameters see also Table 4 of IBruntt et atl ([2004). 
The APASS and VWA methods use the correlation of the 
abundance of neutral Fe lines with excitation potential as an 
indicator of T e ff. The final T e g is found by minimising this cor- 
relation. Hydrogen lines are not considered in APASS and VWA 
methods. It is important to emphasise that in the APASS method 
( IGillon & Magain 2006) the contin uum level is adjusted in each 
iteration. IGillon & Magainl Q2006) also obtained very different 
temperature estimates working with APASS (T e ff =6735+53 K) 
and with equivalent widths (r e jf =6538+44 K). The authors gave 
a favor for the APASS results motivating that the Gaussian 
approximation of the line profiles (used to measure equivalent 
widths) in spectrum of a star rotating with v sin i— 10 km s 
overestimates the equivalent widths. Since we used equivalent 



widths we investigated this possible error source by calculating 
synthetic profiles of 26 Fei lines in the range 6-122 mA with 
our adopted model atmosphere and with Fe abundance of 
log(Fe/A^ t) = -5.05, measured their equivalent widths with 
direct integration, broadened the synthetic spectrum with the in- 
strumental, rotational and macroturbulent velocity profiles, and 
then again measured equivalent widths by Gaussian approxima- 
tion of the convolved spectrum. Inserting both sets of equivalent 
widths in the WIDTH9 code we got an average abundance of 
log(Fe/JV to t) = -5.050 + 0.004 adopting directly integrated 
equivalent widths and of log(Fe/jV tot ) = -5.040 + 0.016 
adopting the Gaussian approximation, as shown in Table |4] 
present in the Online material. No significant dependencies 
indicating any change of microturbulent velocity or effective 
temperature appeared in the second case. Thus we conclude that 
in HD 49933 the use of equivalent widths instead of synthetic 
spectrum should not influence the determination of the model 

parameters. 

The use of both the APASS and the VWA dBruntt et al.1 

2008) methods also leads to a higher T eff for HD 49933 and 
another star HD 32115, for which th e effective temperature 
was derived by IBikmaev et al.l (2002) with the same method 
(Stromgren photometry and hydrogen line profiles) as in the 
present paper. Figure [3] shows individual abundances from 158 
Fe i lines calculated with the model param eters derived in this 
paper (T ef{ = 6500 K, logg = 4.00) and by IBruntt et all (120081) 
(T e ff = 6780 K, logg = 4.24) versus the excitation energy of 
the lower level. For both models the microturbulent velocity, 
derived by us in the usual way, is nearly the same, 1.5kms _1 
and 1 .6 k ms" 1 , respectively. H owever, to be consistent with the 
results by IBruntt et"aTI d2008l) we plot abundances versus ex- 
citation potential for Brun tt's model using v m i c = 1.8 kms -1 de- 
rived in Brun tt et all ([2008). The same dependence was also 
calculated for HD 32115 using equivalent width measure- 
ments from Bi kmaev et al. I (120021) and model a t mosph eres from 
IBikmaev et al.1 (l2002l) ~and from IBruntt et all d2008lh respec- 
tively. As a final check we used the solar flux spectrum. We anal- 
ysed 137 out of 158 Fei lines with the solar model atmosphere 
(Teff = 5777 K, log g - 4.44) also calculated with the LLmodels 
code. There were 21 discarded lines that could not be properly 
fitted in the solar spectrum. Seven lines are blended in the Sun 
which is cooler and more metallic than HD 49933, while the 
rest 14 lines have extended wings and cannot be accurately mea- 
sured. Our analysis resulted in log(Fe/A^ to t) = -4.57 + 0.09 for 
^mic=0.95kms _1 in full agreement wi th the currently adopted 
solar parameters (Aspl und et al.l [2005b . Individual abundances 
derived from the solar Fe i lines versus the excitation energy of 
the lower level are shown on Figure [3] 

In Fig. [3] one immediately notices the negative slope, indi- 
catin g a too high r e ff , adopting the model parameters obtained 
by (IBruntt et al]|2008l) for both stars. This fact together with the 
poor fit of the hydrogen line profiles (see Fig. [U clearly shows 
that both automatic methods based on synthetic spectrum calcu- 
lations, ignoring hydrogen lines, lead to an overestimate of the 
effective temperature by -200 K. With our model we get small 
positive slope that formally corresponds to an underestimate of 
the T e ff by ~30 K. For the Sun zero slope is obtained. 

Our temp erature is close to that derive d with the IR flux 
method by iBlackwell & Lvnas-Gravl (fl998l T eff =6484+45 K, 
and provides a good description of all global observables: pho- 
tometric colors, hydrogen and metallic line profiles. 

Our value for the surface gravity is supported by the ion- 
isation equilibrium for Fe i/Fe n and few other elements, such 
as Sii/Sin and Cai/Can. Even for Ti and Cr, for which 
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Fig. 3. Plots of individual abundances for 158 Fe Lines versus the excitation energy of the lower level for HD 49933 (left panel) and 
for HD 321 15 (right panel) with different model atmospheres. The same dependence for 137 common Fe Lines in solar spectrum is 
shown by crosses in the left panel. 



we have enough lines for an accurate abundance analysis, 
the average abundances from the two ionisation stages agree 
within the error bars. Applying a higher microturbulence of 
1.9 km s -1 , derived from Tin lines, for example, leads to a 
Tii/Tin equilibrium. In the case of chromium about half of 
the analysed Crn lines and one third of Cri lines have the- 
oretically calculated oscillator strengths which may influence 
the final abundance results. The atmospheric parameters de- 
rived for HD 49933 are similar to another solar-type pul- 
sator, Procyon, which has r eff =65 12+49 K, log £=3.96+0.02 , 
and log(Fe/iVtot) = -4.60 + 0.15 dAllende Prieto et alj l2002h . 
The comparison of the observed and the synthetic Mg i line pro- 
file for the 51 72 and 5 183 A Mgi line s for Procyon is shown 
in Fig. |2] (see Ryabchi kova et al.ll2.008l for details concerning 
the observed spectrum of Procyon and the adopted model atmo- 
sphere). 

5.2. Synthetic photometry 

Since there are no available observations of HD 49933 in the vis- 
ible spectral region, calibrated to physical units, we adopted pho- 
tometric observations extracted from the SIMBAD databasfl 
Table [2] summarises the comparison between observed and the- 
oretical color-indexes of different photometric systems. The the- 
oretica l colors were computed using computer codes bv lKuruczl 
(1993a) modified to read and process high-resolution fluxes 
produced by the LLmodels. The reddening, corresponding to 
E(B — V) = 0, was der i ved us ing analytical extinction mod- 
els by Amores & Lepine] (120051) . Besides the final model with 
T e s= 6500 K, logg= 4.0 (Model 1), we present also the syn- 
thetic photometry for three other models: r e ff= 6500 K, logg= 
3.85 (Model 2); a model with the parameters derived from the 
Stromgren photometry (T ef( - 6550 K, logg= 4.2 5, Model 3), 
and a model with the parameters published by iBruntt et"al] 
(120081) (r eff = 6780 K, logg= 4.24, Model 4). It is seen from the 
Table|2]that the Stromgren and UB V indexes are better described 
by Model 1 . The H/? index shows a poor fit for all models, but 
this does not play a critical role in the present study since the 
profiles of hydrogen lines were fitted perfectly with the param- 



2 http : //simbad . u-strasbg . fr/simbad/ 



Table 2. Observed and calculated photometric parameters of 
HD 49933. The values in brackets give the error bars of obser- 
vations. 



Color 


SIMBAD 


t6500g4.0 


t6500g3.85 


t6550g4.25 


t6780g4.24 


index 




(Model 1) 


(Model 2) 


(Model 3) 


(Model 4) 


b-y 


0.270 


0.2754 


0.2727 


0.2724 


0.2399 




(0.002) 










tn\ 


0.127 


0.1251 


0.1255 


0.1330 


0.1355 




(0.004) 










Cl 


0.460 


0.5123 


0.5502 


0.4400 


0.5126 




(0.003) 










HB 


2.662 


2.7272 


2.7291 


2.7279 


2.7501 


B-V 


0.390 


0.3865 


0.3829 


0.3865 


0.3371 


U-B 


-0.070 


-0.0628 


-0.0457 


-0.0942 


-0.0882 


Geneva 












U-B 


1.235 


1.2561 


1.2896 


1.1950 


1.2285 


V-B 


0.491 


0.5125 


0.5166 


0.5123 


0.5700 


B r B 


0.973 


0.9842 


0.9835 


0.9878 


0.9723 


B 2 -B 


1.391 


1.4142 


1.4152 


1.4102 


1.4252 


Vi-B 


1.224 


1.2423 


1.2461 


1.2420 


1.2962 


G-B 


1.586 


1.6135 


1.6185 


1.6146 


1.6847 



eters of Model 1. The color-indexes of the Geneva system are 
generally better described by Model 1,2 and 3 i.e. models with 
lower T e ff. However, the indexes U-B and B\ - B are better 
fitted with a higher temperature (Model 4). Model 1 and 3 show 
a comparable fit, with the difference that Model 1 fits better the 
logg sensitive cl index, but worse the U-B index, while for 
Model 3 it is the opposite. In summary, we find that the majority 
of the photometric indicators and the present accurate spectro- 
scopic study point to the adopted T e g= 6500 K and logg=4.0. 

6. Conclusions 

Based on ten high-resolution, high signal-to-noise spectra taken 
from ESO archive we carried out a precise spectroscopic anal- 
ysis of one of the primary solar-like CoRoT targets HD 49933. 
All ten spectra were averaged to obtain a single spectrum with 
a resolution R * 115000 and signal-to-noise ratio » 500 per 
pixel at A « 5000A. Using this spectrum we revised the funda- 
mental parameters and the atmospheric abundances in a consis- 
tent way empl oying modern ID LTE stellar model atmospheres 
(LLmodels; Shulv ak et al.ll2004h . The derived set of fundamen- 
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tal parameters and abundances provide a good fit to the available 
observables: 1R and multicolor photometry, pressure-sensitive 
magnesium lines, metallic lines and profiles of hydrogen Balmer 
lines. 

It is shown that the implementation of automatic procedures 
for the abundance analysis using only line spectra without hy- 
drogen lines can result in inaccurate parameters and hence abun- 
dances. 
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Wavelength (A) 

Fig. 4. Samples of the first spectrum obtained February 11th, 
the last obtained February 13th and the final co-added spectrum 
(from top to bottom) around the strong Fe n line at ~5018 A. 



-4.9 
Oh -5-2 



-5.4 



Equivalent width (mA) 

Fig. 5. Fei individual line abundances vs. equivalent widths. 
Abundances are derived for the preferred model parameters, 
r e ff=6500K, logg=4.00, 1/^0=1. 5km s _1 . 
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Table 3. Error sources for the abundances of the chemical elements in HD 49933. Column 3 gives the standard deviation <x abn (scatt.) 
of the mean abundance obtained from different spectral lines (internal scatter); a blank means that the number of spectral lines is 
< 2, hence no internal scatter could be estimated. (Note that these values are identical to those given in Table Q]) Columns 4, 5, 
and 6 give the variation in abundance estimated by increasing r e ff by 50 K, logg by 0.15dex, and i/ m ; c by 0.2 km s -1 , respectively. 
Column 7 gives the the mean error calculated applying standard error propagation theory on the uncertainties given in the previous 
columns, i.e., o-\ hn (tot) = o^ bn (scatt.) + o-\ hn (T eS ) + a\ bn (logg) + a 2 ^ (v mic ). For the computation of o-^ bn (tot) of those ions for 
which the internal scatter could not be measured, we have assumed a priori <x a b n (scatt.) = O.lOdex. 
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Table 4. Comparison between the line abundances obtained for a set of Fe i lines calculated from equivalent widths measured with 
direct integration and a Gaussian approximation. 
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84.42 


-5.07 


5165.4100 


73.90 


-5.05 


74.55 


-5.04 


5171.5964 


98.20 


-5.05 


99.20 


-5.03 


5184.2661 


22.80 


-5.05 


23.24 


-5.04 


5191.4550 


98.80 


-5.05 


98.57 


-5.05 


5192.3442 


107.70 


-5.05 


106.98 


-5.06 


5194.9418 


82.90 


-5.05 


84.00 


-5.02 


5602.9451 


67.20 


-5.05 


68.04 


-5.03 


5615.2966 


24.00 


-5.05 


24.30 


-5.04 


5615.6439 


122.60 


-5.06 


120.95 


-5.08 


5618.6327 


15.70 


-5.05 


16.10 


-5.03 


5620.4924 


6.10 


-5.05 


6.24 


-5.04 


5624.0220 


7.80 


-5.04 


7.74 


-5.05 


5624.5422 


72.50 


-5.05 


73.44 


-5.03 


5631.7310 


6.30 


-5.04 


6.39 


-5.04 


5633.9465 


28.50 


-5.05 


29.03 


-5.03 


5638.2621 


31.00 


-5.05 


31.58 


-5.03 


5640.3070 


6.20 


-5.05 


6.33 


-5.04 


5641.4340 


17.50 


-5.05 


17.88 


-5.03 


5655.4930 


32.60 


-5.05 


33.20 


-5.03 


Mean + rms 




-5.050±0.004 




-5.040±0.016 



